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Abstract 

We investigate one-particle semi-inclusive processes in lepton-hadron scattering. In unpolarized 
scattering order Q^ 1 corrections appear only when transverse momenta are detected. We consider 
the twist two and three matrix elements and calculate the semi-inclusive structure functions in 
terms of quark correlation functions. We find that at twist three level not only the standard quark 
distribution and fragmentation function contribute, but also two new transverse "profile functions". 
We demonstrate explicit gauge invariance of the hadronic tensor at the twist three level. The results 
of our approach are used to calculate expressions for some cross sections for semi-inclusive processes. 
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1 Introduction 



In this paper we will investigate unpolarized one-particle semi-inclusive lepton-hadron scattering in 
the deep inelastic scattering ( DIS) limit. We will present the analysis in terms of quark correlation 
functions. These objects encode the non-perturbative parts of the process, which we as yet are not 
able to calculate other than using models. We will present an analysis which is complete up to and 
including 0{Q~ 1 ) where Q 2 is the negative square of the virtual photon momentum. In this order 
we will need to introduce four different projections of the quark correlation functions, called "profile 
functions". The fourier transform of two of these are respectively the quark distribution function and 
the quark fragmentation function which are known from the naive quark parton model at leading order. 
The two new functions arise in 0(Q^ 1 ) and do not have simple parton model interpretations. The 
treatment of (intrinsic) transverse momentum in this process will turn out to be essential for these new 
profile functions. Explicit gauge invariance will be maintained throughout the calculation, as we will 
discuss now in more detail. 

As in all lepton scattering experiments in the one-photon-exchange approximation the semi-inclusive 
cross section will be essentially given by a contraction between a leptonic tensor and a hadronic tensor. 
All the interesting physics resides in the hadronic tensor. For the one-particle semi-inclusive cross 
section the hadronic tensor is given by (see section 2 for exact definitions and conventions) 

= 7A4 / ( 2 t)32E x I d4xel9X < P \Mx)\Px, Ph ><Ph,Px\MO)\P> ■ (1) 

Here q is the virtual photon momentum, P and ph are nucleon and observed hadron momentum and 
Px is shorthand for the ensemble of unobserved hadrons in the final state. Also the phase space integral 
over Px should be understood to be a summation over all n-particle states of all types. We denote this 
tensor by the diagram in figure [l|a. In terms of quarks and gluons we know the photon current, which 
only couples to quarks. This gives us an expression for the hadronic tensor in terms of quark fields 

2MW U „ = ^ / J*?* I d 4 xe^ x 



^ (2vr) 4 7 (2^) 3 2£ x 

x <P\^(x)^(x)\P x , Ph ><p h ,PxmOh^(0)\P>, (2) 

depicted by the diagram in figure |l|b. Flavour degrees of freedom are suppressed throughout. Because 
of the colour and electromagnetic gauge invariance of these quark currents the hadronic tensor is 
gauge invariant and q^W^ = 0. We will make a diagrammatic expansion of which the first term is 
the "Born" diagram, i.e. the naive parton model, and the following diagrams have successively more 
gluons or quarks participating. After a specific choice of gauge this diagrammatic expansion will be 
the same as an expansion in Q^ 1 . The relationship between diagrams with successively more gluons 
or quarks participating and their order in Q^ 1 or "twist" is discussed at great length by Jaffe and Ji in 
The separation of a diagram in a hard scattering part and soft correlation functions is discussed by 
Ellis, Furmanski and Petronzio (EFP) [|2| and by Collins et al. in |||. Qiu [Q] and Ji || among others 
discuss this approach for inclusive scatterings. Several authors have used these methods for polarized 
lepton-hadron and polarized Drell-Yan scattering [Q, §]. In the same spirit we will use this method for 
semi-inclusive lepton-hadron scattering. The Born term, shown in figure |2|, is given by 



2MW B ^ LU = — \t [ d A xe 

(27TJ 4 J 



x 



< P\ ^ (x)^ (0) \P X 0|Vj {x)a{a h ^ k (0) |0 > 



+ <P|^(x)^ fc (0)|P><0|^(x)aT ah ^(0)|0> • (3) 
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Several comments are in order here. We now have two separate, unobserved final states which are 
integrated over using completeness, represented by the top and bottom parts of the diagram in figure 
^. In the top part of the diagram the additional hadron in the final state appears. In appendix A 
it is shown how to handle this when integrating over the final state. The result is that we are left 
with a hadron number operator. In the diagram we see the separation in quark correlation functions 
representing the soft part and the two gamma matrices that denote the hard photon-quark scattering. 
Furthermore we note that eq. (||) contains contributions of quarks and antiquarks. In the rest of this 
paper we will not explicitly account for the antiquarks since they behave exactly like the quarks in a 
distinct sector. Of course the numerical values of quark and antiquark correlation functions can be 
quite different. Finally we note that we will be concerned with hadrons produced from the struck quark, 
i.e. belonging to the (isolated) current jet. This will pose restrictions on the scaling variables x B and 
z, as will become clear later in the discussion of kinematics and factorization. 
Evaluating eq. (|^) we will find the following contributions 

5 + (i 



Wj(Q)=»S, +< 1 -^+O( 7 ^). (4) 



where q^w^® = but q^w^ / so we explicitly loose gauge invariance in 0(Q^ 1 ). This is not 
surprising since the quark fields in the matrix elements in eq. (|3|) are taken at different space-time 
points, so this matrix element will in general not be invariant under a local gauge transformation. 

We will discuss how one can restore gauge invariance in such a situation before returning to the 
diagrammatic expansion. Gauge invariance of a bilocal matrix element ^(x)ip(x') requires a link 
operator ip(x)L(x,x')ip(x'), which can be written in the form 

L(x, x') = V exp(ig J d 4 y £ Q (x -y,x'- y)A a (y)), (5) 

where A is the appropriate gauge field Q. Gauge invariance is restored when ^ satisfies 

9 -C(x -y,x'-y) = 5\x' - y )-5\x-y). (6) 



dy< 

The path ordering V is then defined along the flow lines of the "velocity" -field £ a . A special choice of 
£ leads to the well known line integral 

L(x,x>) =Vexp(ig (* dy a A a {y)). (7) 

Jx' 

So if one expands in orders of g, which is in our case also an expansion^ in orders of Q^ 1 one concludes 
that restoring gauge invariance requires addition of a term 

i$(x) J A - y, -v)Mv)*l>(P) (8) 

to the Born term. 

We return to our diagrammatic expansion, and inspired by eq. @, consider in addition to the Born 
term the four diagrams of figure [3|. They have one fermion propagator in the hard scattering part and 
their leading contribution is 0{Q~ 1 ). Generically 

W$(Q)=w*f-± + o(±). (9) 



Q \Q 

The essential part of our proof of gauge invariance will be as follows. We will show that in these 
diagrams the soft parts are matrix elements of operators of the form in eq. (||) and in this way identify 

This is only true in a specific gauge, we will discuss this later on. 
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£ a in our case. Then we will show that this £ Q respects the condition eq. (|6|) up to 0(Q~ 2 ). We claim 
that all other diagrams with extra gluons either give rise to the logs of the radiative corrections, or are 
suppressed by powers of Q, or ultimately vanish in a suitably chosen gauge. 
Summarizing we find 

w M „ = <° + « + • 1 + o (ij ) (io) 

and q^W^y = 0(Q~ 2 ) in our calculation. After explicit evaluation of these five diagrams we express 
W^v in four unknown profile functions, fourier transforms of specific projections of the quark correlation 
functions. 

The rest of the paper is structured as follows. In section 2 we discuss the generalities of the semi- 
inclusive cross section, kinematics, structure functions and inclusive reductions. In section 3 we will 
first look at the kinematics of the relevant diagrams and discuss what kind of assumptions are made 
related to the factorization into hard scattering and correlation functions and make a twist-analysis, 
then we will define the basics objects of which our diagrams are composed and show how they relate to 
the four profile functions which characterize the process. Then in section 4 we will present the actual 
calculation of the hadronic tensor and demonstrate gauge invariance. Finally in section 5 we present 
our results as expressions for various cross section using these four scalar functions. 



2 Semi-inclusive deep inelastic lepton hadron scattering 

We will consider semi-inclusive processes of the type eH — > e'hX, where H is a hadronic target with 
mass M, h is a hadron with mass my t detected in coincidence with the scattered electron and X is the 
rest of the final state. The momenta are defined in figure [l], the angles in the target rest frame (TRF) 
are defined in figure 0. In a one-hadron semi-inclusive process we can form four independent invariants 
q 2 , P ■ q, P - ph and ph ■ q- This includes the invariants that one is familiar with in inclusive scattering: 

Q 2 = V^^Wg), (11) 

v = ^— — ^ T = F E — E 1 , (12) 

and the ratios 

x v-^~- (13) 

Xb ~2Mu' y -p-k~E- {U) 

For semi-inclusive scattering one or more particles (=momenta) in the final state are measured. The 

additional invariants are fixed by the energy E^ and the component of ph along q, indicated by Ph\\, 

or equivalently by Eh and p|j_. In the TRF they do not depend on the azimuthal angle 4>h- Note 

that if more than one momentum in the final state is measured the relative azimuthal angle appears 

in the invariant p\ ■ p2- In particular in the TRF one has p li _ ■ p 2 ± = \Pi±\\P2±\ cos(0i — <fo)- l n DIS 

processes the ratio 

P -Ph trf E h 

z = ~b = — ( 14 ) 

P ■ q v 

will be useful. We will consider in this paper the limit of deep inelastic scattering, in which Q 2 — > oo 
and v — > oo, keeping the ratio xb fixed. In this limitf] (defining z along — q), q~ — ► oo. One then 
immediately sees that for a particle produced in the current jet the other invariant involving ph and q 
with ph on mass shell is: 

2Ph-q _ Ph 1 + _ Ph 



q 2 q q + q 



(15) 



v-2 
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The cross section for the semi-inclusive process is given by 



d 3 p h dQdE' Q 4 E 

or 



W (16) 



2E h da ^My L ^ ^ 



d 3 p h dx B dy Q A 
where L^ v is the spin averaged lepton tensor 

L[iv = 2k^k u + 2k I/ k fM — Q 2 gp, v (18) 

and the hadronic tensor is the product of non-elementary hadronic currents, 

2MH V = 7A4 / \ 2 d Jf2E x < P \M0)\Px,Ph><Ph,Px\M0)\P> 

x (2n)H A {P + q-p h -P x ). (19) 

The integral over Px indicates a complete summation over all multiparticle states of all hadron types. 
We have split off the invariant phase space for the detected hadron, which can also be written: 

f£ = = \^izd P l^ h . (20) 

ZEh *Ph\\ ^Ph\\ 

The most general gauge invariant, Lorentz covariant hadronic tensor for the electromagnetic process 
(parity and time reversal invariant) can be decomposed in four structure functions Wi(x,Q 2 , z, p|j_) 
I: 



. Ph±tiT u + T^p h±u Ph±/iPh±u , 0I , 

+ W 3 + — ^— W 4 , (21) 



where 

T, = P M ~ (22) 

p£± T =(0, Ph±,0). (23) 

As in the inclusive case, it is convenient to introduce photon polarizations. The structure functions can 
be expressed in terms of: 

q 2 

W L = e L ■ W ■ e L = -Wi + ^ W 2 , 

Q 2 

W T = \ (e x ■W-e x + eyW-e y )=W 1 + -^ W 4j 

Wlt cos(^) = -(e x .W-e L + e L -W-e x ) = -2^^^ cos(^) W 3 , 

Wtt cos(20 h ) = \ (e x ■ W ■ e x - e y ■ W • e y ) = cos(2fo) W A . (24) 
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The e's here are the standard transversal and longitudinal photon polarization vectors. The tensors 
multiplying these structure functions in 



1 



WUP,Ph, Q) = PuuWl + P^Wt - ^PuJWlt + P'J Wtt 



2 ^ 



[IV 



are 



P T 

[11/ 

oTT 



T^T V 

J<2 ■ 



y " r g2 " r p2 ' 
Ph±yT u + T^p h ± v 
T \ Ph±\ 

^ _ Q^Qv _ PjiPy r, Ph±nPh±u 



y2 



PL 



They can be used to project out the structure functions W, 

>V LT 



P^W^, 

P^JW^, 
2 



as they are orthogonal to each other. 

In DIS it is convenient to use the combinations Hi = Hi(x,Q 2 ,z, p^±) 



2m\ 



2zHi = M(Wi + ^W 4 ) = MWt, 

2zH 3 
2zH A 



2mf i q 2 



-^-(W l + Wt), 



2 

Mx B -^W4 
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For later purposes it is also useful to introduce 

2zH L = MW L 



2z 



-Hi +11 + —z- 



Q 2 \ n 2 



v l ) 2x, 



DIS 



2z 



-Hi + 



H2 

2x F 



(25) 



(26) 



(27) 



(28) 



(29) 



The cross sections are easily obtained from the contraction of leptonic and hadronic tensors. Without 
any approximations one obtains 



2E h da 
d 3 p h dfldE' 



da 



Q 2 1 



dQJ M q 2 e 
+e Wtt cos(20 h ) + 



W T + eW L 



le(e+l) 



Wlt cos((/>h) 



(30) 
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where 



and 



da 



2 TP/2 



M 



Aa z E 

,2 



cos l- 



1 + 2^ tair ( - 



(31) 



(32) 



It is also possible to observe a second particle in the final state. One then finds by a similar analysis 
two extra structure functions associated with the extra particle, W 2 T and Wf-y, which enter the cross 
section in the same way as the one particle structure functions, (see [||]). 

Up till here all expressions are completely general. Now we want to specialize to DIS events. Our 
attention will be the current jet, considering particles produced with bounded transverse momenta, i.e. 



™l+ PhX^PhW - 



El 



(33) 



The cross section may be written as (still exact): 



da 



4vr a 2 ME E h 



dx B dzdyd p^dtph 



Q 4 
Ph±\ 



Q 



Ph 

(2-//) 



x B y H\+[l-y 



Mx B y 
2E 



1-2/ 



Mx B y 



+ 



y 



Mx B y 
2E 



Q 2 



i + 



2Mx 1 
yE 



yE 

cos{2(j) h )TiA 



(34) 



which becomes in the DIS limit 

da 



Dis kna 2 ME 



dx B dzdydp h ±d(ph 



+■ 



Q 4 

Ph±\ 
Q 

,2 



x B y 2 Hi + (l-y)Hi 



(2-y)yJ{l-y)cos(<t> h )H 3 



cos(2^)W 4 



(35) 



Integrating over the azimuthal angle only leaves the first two terms. Often one is only interested in the 
z-behaviour of the structure functions, so defining 



Hi(x B ,Q 2 ,z) = - Jdp 2 h± dcf> h Hi(x B ,Q 2 ,z, p 2 h± ) 
= 7T Jdp 2 h± Hi(x B ,Q 2 ,z, p 2 h± ) 



we are left with 



z s y 2 #i + (1 - y)H 2 



(36) 
(37) 



da dis 8vr a 2 ME 
dx B dydz Q 4 

In the analysis of experimental data, factorization of the structure functions is taken to mean that 
one can write the structure function H 2 as a flavour sum over products of quark distribution functions 
and quark fragmentation functions (K 



H 2 {X B , Q 2 , Z) = H 2 (x B , Z) = e l X B fq/H{x B )D h / q {z) + fq/H(x B )D h/ g(z) 



(38) 
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where q indicates quarks of different flavours. In this paper we show the validity of this factorization 
in leading order in Q, and its extension to 0{Q^ 1 ) and secondly we derive an expression for the 
fragmentation function in terms of quark correlation functions. Furthermore one obtains 



2x B Hi(x B ,z, pl ± ) =n 2 (x B ,z, pl ± ), 



(39) 



as H L -> in the DIS limit. 

Finally some remarks on the connection to inclusive scattering must be made. The hadronic tensor 
for inclusive scattering W /Jtl/ (P, q) is given by (see appendix A): 



(n h (P,q))W^(P,q)= J ^W^(P,q,p h ), 



(40) 



where (rih(P,q)} is the average number of particles produced of type h for given P and q. Combining 
the definitions we then have 

(n h (x B ))F 2 (x B ) = J dzH 2 (x B ,z), (41) 
where F-i are the standard dimensionless scaling structure functions in inclusive IH . 



3 Factorization and correlation functions 
3.1 Factorization and twist analysis 

In this section we will discuss the diagrammatic expansion. We start by considering the kinematics. It 
is convenient to use lightcone coordinates for this. With the four vector notation p = [p~,p + ,p±] with 
p ± = (p° ±p 3 )/\/2 we denote the momenta of the target hadron (P), virtual photon (q) and final 
state hadron (ph) as 



P = 



x B M 2 A 



Ay/2 ' x B V2 



,0, 



Q 2 A 

:> 7^' U -L 



Ph 



Ay/2' y/T 
zQ 2 M p 2 h± + ml) 
Ay/2' zQ 2 V2 



Ph± 



(42) 



Here A is a free parameter fixing the frame among those that are connected by boosts along the 
photon-quark direction. For instance A = x B M corresponds to the target rest frame (TRF), while for 
p h± = 0, the choice A = (zQ 2 )/mh corresponds to the outgoing hadron rest frame. A convenient 
infinite momentum frame (IMF) is found by choosing A = Q. Using the momenta P and q, and the 
momentum T = P — (P ■ q/q 2 ) q, one can construct the momenta 
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Q 



A 



Q 



Ay/2' QV2 

Q A 
Ay/2' Qx/2 



,0± 



,0± 



(43) 
(44) 
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and the two lightcone unit vectors 



n% = ML-J (Z* - (45) 

fit = (-^\(f» + n (46) 

with h + • n_ = 1. Note that in our results the invariants x B and z are used as the ratios of lightcone 
coordinates. Thus the replacements 



q + 2x 



D 



P + / 4M 2 x 2 



(47) 



would provide the necessary socalled hadron mass corrections. Furthermore we assume that for the 
current jet events we consider, no transverse momenta growing like Q appear. We consistently neglect 
terms of order transverse momentum divided by Q 2 . Also x B and z should not be too close to zero. 
To be precise, x B Q and zQ must be at least larger than some hadronic scale. 

Diagrams contributing to the process will be separated into soft hadronic matrix elements and a 
hard scattering part. The forward matrix elements involving the external particles (momenta P, ph) 
are given by untruncated (except for external lines) Green functions. Given a renormalization scale 
for these Green functions, they are calculated as 2-particle irreducible in the quark anti-quark channel. 
The evolution between different renormalization scales is then provided by the logarithmic corrections 
discussed below. It is assumed that for all of the Green functions the parton virtualities are restricted to 
some hadronic scale, i.e. they vanish with some inverse power of the virtualities. The relevant diagrams 
that are considered in this paper have been given in figures ^ and ^. There are other diagrams with 
one gluon such as those of the type in figure || with the photon-quark-quark vertex and the gluon- 
quark-quark vertex interchanged. Those can be absorbed in the soft parts and will not be considered 
explicitly in our analysis. Diagrams with one gluon which are important are those shown in figure 
||. They contain in the hard part partons in the final state and will produce logarithmic corrections 
from the collinear divergencies^. This leads to the evolution of the profile functions f + ,D~, also 
known as the Gribov-Lipatov-Altarelli-Parisi equations, as has been proven explicitly by Ellis, Georgi, 
Machacek, Politzer and Ross [11]. The non-logarithmic terms in their calculation are the gluon jet 



events which appear in order a s and which have been calculated in a naive parton model by Konig 
and Kroll in ref. [12|. We expect that from diagrams like those in figure which have in addition 



gluons in the final state one will obtain Q~ l logQ 2 contributions. In this paper however, we will not 
discuss the logarithmic contributions any further. We do note that in the leading order we can calculate 
the logarithmic corrections in our approach and that the evolution equations are not altered by the 
observation of transverse momenta. Finally there are diagrams which have a gluon across between the 
hadronic parts. With the above assumption on the parton virtualities one immediately sees that they 
will vanish as some power of Q~ 2 . The general behaviour of the contributions of diagrams with any 
number of partons emerging from the hadronic matrix elements is discussed below. 

It is illustrative to first consider the kinematics of the diagrams in figures ||| and fj. We assign the 

3 The quark self-energy and the photon-quark-quark vertex correction provide singular terms needed for the mass 
factorization. 
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four momenta 



V 



Pi 



t{p 2 +pl) J_ 

(Q 2 A{k 2 + kj) 
Ay/2' (Q 2 V2 

(p?+P?jl) (1- 



P± 



(l-^)Ay/2' V2 
(l-Ci)Q 2 Af&f + fcfJ 



Ay/2 '(l-Cl)Q 2 ^' 



1J_ 



(48) 



Now we apply four momentum conservation p + q 



is equivalent to jf 1 



x B P + and fc 



to derive £ = £ = 1 in the limit Q 2 — ► oo, which 
= zpfr- For later use we also give expressions for 



the momenta in the quark propagators in figure 

Pl+q 



(pi + q) 2 

h - q 
(fci - q) 2 



-l 



.4 



= = -Pi± 

1 -vl -feii. 

aTI 7 CiQ 2 V2' W 



(49) 



It is important to note the suppression of the transverse components relative to the other components. 

The objects that enter the diagrams in figures || and fj are depicted in figure [|[ They are partly 
given by untruncated (except for external particle lines) Green functions for the blobs and partly by 
ordinary QCD Feynman rules. Suppressing flavour and colour indices we have in the hard part the 
following rules, 



5a 
5b 

5c 



i 

f 



while the correlation functions are given by 

1 



5d 


T iM = 2 


5e 


Dji(k) = - 


5f 


F?j(p,Pi) 


5g 


M«(kM) 



d^x e~ ip - x <P\iP j {xW i {0)\P>, 



1 



4z(2tt) 4 
1 



d"x e lk ' x <0\ip j {x)ala h ip i (0)\0>, 
d 4 xd 4 y e - ip - x e i( P-Pi>y <P\^ j (x)A a (y)i; i (0)\P>, 



2(2tt) 8 
1 



(50) 
(51) 

(52) 



(53) 
(54) 
(55) 



4z(2ttY 



d*xd*y e m ' x e 



ik-x —i{k—k\)-y 



< 01^)^)4^(0)10 > • (56) 



We note that time ordering is not relevant. Using the QCD equations of motion |13[] all correlation 
functions can ultimately be related to functions containing ip + = (7 _ 7 + /2)ip and A 1 - , which have 



c-number (anti)-commutators at equal lightcone time x + [14]. 

We will now discuss the Q-dependence and give the general classification of all possible diagrams 
with more than one parton emerging from the hadronic matrix element according to the powers of Q. 
We use an extension of the method developed by EFP for inclusive scattering. They have proven that 
for inclusive scattering this method is equivalent to the Operator Product Expansion (OPE). Because 
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of the hadron in the final state we do not have an OPE for semi-inclusive processes. We can, however, 
follow the method of EFP for semi-inclusive processes. 

Consider an observable, say one of the structure functions Hi in the expansion of W^„. Separating 
the hard scattering from the soft pieces we have schematically 

0(P,q,p h ) = J2°n,m(P,<l,Ph), (57) 

n,m 

where 

„ n m 

O n ,m(P,q,Ph) = / nil dpidkj f(P,pi)T(q,pi,ki)D(p h ,ki). (58) 

i j 

Omitting the diagrams with partons in the final state, a general term O n/m , with n + m partons 
participating in the hard scattering is depicted in figure [/]. To investigate the behaviour of the diagrams 
for high Q 2 it is convenient to work in the particular infinite momentum frame A = Q. Under the 
assumption that all parton virtualities must be restricted to some hadronic scale we can show as seen 
explicitly for the momenta appearing in figures § and ||, that kj and pf ~ Q while k~j~ and p~ ~ Q~ x . 
The transverse momenta are all of order unity. For the truncated piece T one can convince oneself 
that in the limit Q — > oo, the Q-dependence is determined by dimensional arguments. Assume that 
the canonical dimensions of the (in general nonlocal) product of operators appearing in T and D are 
dx and do- The dominant hard scattering part multiplying T and D then must behave as 

r( Ph k j ,q)~Q c - d T- d °, (59) 

where c is a constant depending on the specific observable. The behaviour of T is determined by 
its Lorentz structure and the specific component that one considers. As in the matrix element every 
Lorentz-index contributes as P M one finds after collecting the powers of Q in the frame A = Q where 

P + ~ Q, P~ ~ M 2 /Q and P 1 - ~ M 

f ~ M dT - hT Q hT . (60) 

After determining hx in this way, the twist of the (nonlocal) operator in T is defined as tx = dx — hr- 
Similarly, the behaviour of D, where p^ ~ Q, p\ ~ rn^/Q and p h± ~ m^, is given by 

D ~ m d h D ~ hD Q hD (61) 

and the twist is defined as to = dp — hp. Note that we (following Jaffe) discuss twist for an operator 
which is in general nonlocal. Moreover, different components (defining different profile functions) will 
have different twist. For example V ; (^)7 + ^/ ; (0) will have t = 2, but ip(x)j~ip(0) has t = 4. The twist 
for the + + ...+ component of a symmetric traceless operator is identical to the 'twist' as used in 
the OPE approach for inclusive DIS. The name twist is still used in our context as the above leads for 
O n ^ m to the behaviour 

/ 1 \ t T +t D -c 

0{P,q,p h )~^-) . (62) 

The leading Born diagram starts off with bilocal twist two operators in T and D, but also contains 
twist three and twist four pieces. Diagrams with, for instance, an extra pair of quark fields in T or 
D are suppressed by a factor Q~ 2 as tx or to is increased by 2. For gluons the situation is more 
complicated, since addition of a gluon does not change the twist. By choosing the gauge A + = for 
the lower and A~ = for the upper blob we force the gluon contribution to add at least one unit of 
twist. Let it be clear however that one can add in general an arbitrary number of longitudinally gluons 
to every diagram. Only after a suitable choice of gauge these contributions will vanish. 
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3.2 Correlation functions and profile functions 



The profile functions will be defined as fourier transforms of projections of T and D. Explicit definitions 
of the profile functions are given in appendix B. With these definitions we find 



dp- Tr [T(p)L 



-X B 



p+ 



q 2Mx B 

-Q-Q-<*»>P±)> 



(63) 



dp- Tr [T{p)^] p+=XBP+ 



q_ 
Q 



2x B p p L 



+ 



Q 



-f ± (x B , P± )+0 



dk+ Tr[D(k)] k - =p - /z 
dk+Jr[^D(k)] k . =p - /z 



q + 2m h 



(64) 
(65) 



Q 



(f»+r) d-(z, p h± -zk ± ) 



+ 



+ 



(66) 



where / + and D~ are twist two profile functions and e, f- 1 , E, D 1 - are twist three profile functions. 
These are all scalar functions. The functions / + and e are the same as defined by Jaffe and Ji in [0]. 
The function D~ is the same as defined by Collins and Soper in [15|. The functions / and D 1 - are 
new. 

In order to handle the three-parton correlations M and F we define the functions M. and T, 
completely analogous to M and F but with the replacement A a — > D a = d a — igA a . With partial 
integration we proof the relations 



J d*pi Tij(p,pi) = -ip a Tij(p) -ig J d 4 px F^(p,pt), 



d 4 h MfAk,h] 



-ik a Dij(k) -ig / d 4 h M^Jk,h). 



(67) 



As argued before, the three-parton correlation functions F and M or T and M. will be twist 3 in 
a suitably chosen gauge. We will now show that the associated leading profile function is in fact 
expressible in the twist three profile functions appearing in eq. ( |64D and eq. (|6^). For this we use 
the QCD equations of motion g |[ [TJ]. We define ip± = P ± ip where P ± = 7 =F 7 ± /2. It is then 
straightforward to show that ip + and A 1 - are dynamically independent fields in the gauge A + = 
1 17 ] . This means that parts of the three-parton correlation functions corresponding to ^_ and A~ can 
be written as correlation functions involving more than three partons and are thereby guaranteed to 
contribute at twist > 3. From this simple argument we know beforehand, that if we are interested in 
twist 3 contributions (and we are!) then we only have to consider Trpy + F a (p,pi)], with a = 1,2. 
For the M functions we find twist 3 for Tr[y ~ M a (k,k\)] in the gauge A~ = 0. One draws the 
same conclusion from the twist analysis of the previous section. As = \^2ifj\_ip + we conclude 

that the leading contribution (in twist 2) comes from ip\.A + t/) + . When choosing A + = the leading 
contribution is then at twist 3 and given by ijr, A^ij)+. 

In the A ± = gauge we deduce the equations of motion 



D a {x) 1 ± ip{x) = d ± 1 a ^{x) 



(68) 
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This is only true when applied between unpolarized nucleon states. Applying these equations we find 
for the relevant transverse components (a = 1,2) 



Az(2tt) 



dk+d A x e ik '- x Jr <0\^(x)iD a (0)ala h ip(0)-/~\0> 



(69) 

^£>\ (70) 



Here k' is the quark momentum k obtained after a Lorentz transformation to the frame where p h± = 0, 
k'j_ = k± — Ph±/z as explained in appendix B. We have to transform back to the original frame with 
Ph± 7^ 0- Note that because of this transformation 



4z(2vr) 



dk + d A x e ik ' x Jr < Q\ij}{x)iD- (0)0^0^(0)7" |0> 



k~D~ 



(71) 



will contribute to the transverse three-parton piece. This gives us as the result in leading order which 
we will use in the next section 



dp- I d A Pl Tr[iTYV = _^ 
dk + I d 4 h Tr [iM a ^] k - =q - 



Q 



(t^ - r)pt x B r + v[q2 



1 



Q 



(T M + qt- 



z z 



+ 



(72) 
73) 



4 Calculation of the hadronic tensor and gauge invariance 

We present the calculation of the hadronic tensor W from the diagrams of figures [2] and [||. It consists 
of three parts, W = Wb + Wf + Wm, where Wb is given by the Born diagram of figure ||, Wf by 
the one-gluon diagrams FJ|a and [|b and Wm by figures [|c and [3|d. 



Az 
M 
-\z 



d 4 pd 4 k 5 4 (p + q- k)Jr [^T(p)j u D(k)] 



(74) 



d A pd A k8\p + q-k) J d 4 Pl (Tr[ 7 ° ^ gF a (p, Pl )^ D(k)] 

J W + §2 l a D{k)^gF a (p! , p)] ) , (75) 



-Az 



+Tr[ 7 Q 



d [ P d l k ^(p + u-k) jd% (Tr[ 7 M ^ 2 7 a T{p)YgM a (k, h)] 



(h ~ q) 



rfgMaikiM-fTfr)]) 



(76) 



It is not necessary to consider the colour charge operators at this point, as they will disappear again when 
the gluon fields are eliminated using the QCD equations of motion. As discussed in the introduction 
we will show that Wf and Wm generate the first order of an expansion in g of the linking exponential, 
needed to render the Born term gauge invariant. We will discuss only the F part as the M part is 
completely analogous. The required term for gauge invariance is 



KV Link,F 



Az 
M 



d A pd A k 5 4 (p + q-k)Jr 



Ikl 



2{2ir) 4 



d A x e~ ip - x 



x / d 4 y <PW i (x)$ lk , f (x-y,-y)A a (y)MQ)\P> % 



Az(2n) 



d A z e ik - z <0|Vy(^Ka/^ fc '(0)|0> 



(77) 
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where we have the condition on £ Q 

d 



dy 

By comparing Wp U with W^ k F we identify 



£k>k,jj>(. x -yi-y) = s jf 



d 4 pi 

(27T) 4 



7 



fa + 



Ap-pi)-v 



k'k 



d 4 p\ 
(2vr) 4 \i> x + i 



7 



(78) 



J(p-pi)-(x-y) 



33' 



We calculate the derivative of £ Q and use the identities 

1 



(p-pi) c 



fa + i 
{p-pi) a ii 



la 



fa + 



fa + i 



fa + 



1, 



(79) 

(80) 
(81) 



which are essentially Ward-Takahashi identities for the hard scattering vertex. In leading order we use 
the kinematics of eqs. (|8|) and (fw|) to find 



^k' k , 3 A x - y> -y) = 6 Jf p kV(y) - ^ k p- r 5 4 ( x - y) + o ( 



(82) 



Since for the D correlation function the leading parts are given by V- = P~ip fields in the matrix 
elements and (P~) 2 = P~ , we see that in leading order our £ a fulfills condition (|7£|). From this and 
the analogous reasoning for the other two diagrams we conclude that the four diagrams of figure 4 are 
enough to render our calculation gauge invariant up to 0{Q~ 2 ). Note that both in our derivation of 
the linking operator in the introduction as in our determination of the leading diagrams we have not 
considered the contributions of multiple longitudinal gluons. However, we have proven gauge invariance 
for the 0{g) term which in the gauge A + = 0, A~ = respectively will be the leading term. In our 
calculations we will also not need to consider these diagrams as we will choose these gauges. 

Now that we have convinced ourselves that our starting expression is gauge invariant we return 
to the calculation of the diagrams. We will expand the correlation functions in the standard basis of 
Dirac matrices. Because of time reversal and parity invariance together with hermiticity only 1 and 7 M 
can contribute in unpolarized electromagnetic scattering. Furthermore we can proof with time reversal 
invariance and hermiticity that Tr [7 M F a (p,pi)]=Tr \y^F a (pi,p)\ and the same for M a (k,k\). The 
tensors reduce to 



4z_ 

M 
-4z 



J d 4 pd 4 k 5\p + q - k)±Tr [ 7(T T] Tr [ lx D] Sjf 
d 4 pd 4 k 8 4 {p + q-k) 



M 



-Az 

w 

d 4 h 



+ 



(83) 



(84) 



d A pd A k 5 4 {p + q-k) 

(fci - q) f 

(h - qf 



[ ^Tr [ la T] Tr [ lx M a {k, k x )\ {s^ x + S a ' wX ^} ) , (85) 



where 



qUcfvX 



Tr 
Tr 



ryr ryt~ ry ry ry 



(86) 
(87) 
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Contributions of the Dirac unit matrix either give rise to an odd number of gamma matrices or are 
suppressed by Q^ 1 . Since the fermion propagator in the frame A = Q goes as Q^ 1 we only have to 
pick out the order 1 contribution of the rest of the integrand. The leading contribution is given by 

for Wp and 

for . Because 7+7+ = 7^7^ = and { , y ± ,j ± } = this means that p = — in the first case and 
p = + in the second. This in turn implies that the fermion propagators are in fact constants in leading 
order as can be read off from eq. fl49| ) in the frame A = Q 

= iff-gY\ (90) 



(pa + g) 2 2Q V 

Now we can take the propagators outside the fci-integration and use the relations (|67|). Furthermore 
we decompose 5 A {p + q — k) as 5{p + + q + ) 5(q~ — k~) 5 2 (p 1 _ — k±). We finally find 



z 
AM 



J dp + dk-d 2 p ± ^Tr [ la T] Tr [ 7x D] S% av \ (92) 
(Pi + g) 

(pi + q) 

Tr [ 7ct T] Tr [ lx D] P±a ) [ S a ^ x + ) j , 



W? = ^Jd P Uk-d 2 P± I^U/^VTr^] Tr[ 7A ,D] 



(93) 



W A T = -^Jd P + dk-d 2 p ± (Tr[T 7g ] Jd%Tr[iM alx ] 



(h - q) 2 

- Tr [ 7(T T] Tr [ 7A £>] A; ±Q ) ( + go^A^ j | _ ^ 

This is the contribution of the relevant diagrams, expressed in quantities which we have given as tensors 
in next-to- leading order on the basis {T, q, p±,kj_,ph±} in the previous section. To be precise, we will 
substitute the relations (0), ©, @, © and (f|). Note that f 2 = 1, q 2 = -1, f • q = 

and that they do not have transverse components. The calculation is then very well suited for symbolic 
manipulation with a computer. We have used FORM |1£] to evaluate the total hadronic tensor as 

W"'(P,q,p„,p ± ) = (-«" - 9"") (^Z+O- 



+ [t"t-) (^rD- 

+ (f A, + (^/ + °- - ^/+^) +o (£) , P5) 

where 

)T = / d 2 p ± W^. (96) 
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One can convince oneself that q^W' 11 ' = 0(Q~ 2 ). The hadronic tensor is in fact the hadronic 
tensor for the semi-inclusive process where one observes in addition to a hadron h the transverse 
momentum of the jet, i.e. the transverse momentum of the ejected quark in our approach 



j (e+H->e'+h+jet) a 2 E' 1 ~ 

= __ —L 

d 3 Ph d 2 P± dndE' Q A E 2E h ^ 



(97) 



As mentioned in section 2 this gives then rise to two extra structure functions. We project out the 
leading structure functions (all other structure functions are of order Q~ 2 ), 



W T (.X B ,Z, P h± ,P ± ) 

Wl T (x B ,z, p h ±,p±) 
W 2 LT (x B ,z, p hX ,P±) 



z 1 
— f + (x B ,p ± )D~(z, p hl _ - zp ± ) + C(g2)> 

\Ph±\ 4 r +( 

x(^D ± (z, p h± - zp ± ) - D~(z, p h± - zp ± )) 
\P±\ 4 



z{f + {x B ,P±) -x B f (x B ,p x ))D (z, p h _L-zp ± ) 



Q M 

- f + (x B ,P±)D ± (z, p h± - zp ± 



(98) 



5 Results 

5.1 The quark distribution and fragmentation function 

In this last part we will present the cross sections of several processes in lepton-hadron scattering. 
These cross sections will be expressed in the previously derived quark profile functions. Firstly we will 
examine the properties of these functions. 

The profile function /+ is just the well known quark distribution function, see [19]. The function 



fq/H{XB,P±) = f + (XB,P±) 



2(2vr) 



^— / dx~ d 2 x ± e< q+x +P±- X ^ <P|i/>(x)7+^(0)|P> 
lirr J 



x+=0 



(99) 



is interpreted as the probability of finding a quark with lightcone momentum fraction x B = p + / P + = 
—q + /P + and transverse momentum p^ in a target with no transverse momentum. This can be seen 
when one quantizes the quark fields on the lightcone with x + = (see (lj). As the fields are 
the essentially free fields one can then substitute a free field expansion. The operator in the matrix 
element then essentially reduces to a quark number operator, counting quarks with momentum fraction 
x B and transverse momentum p ± . The integral over x B and p ± comes out as 1, so / + is a probability 
distribution. Upon integration over the transverse momentum one finds the the quark distribution 



fq/H( 



1 

4tt 



dx 



iq • x 



<P\^(x)~f + ip(0)\P> 



x+=X 



(100) 



With this form one can easily prove sum rules expressing probability and momentum conservation 
at the quark level. As a check on the answer one can consider a quark target. One expects that 
the probability of finding a quark in a quark reduces to a delta function in x B and p ± . An explicit 
calculation, substituting a free quark state for the target state \P> leads to 



f g / q (x B ,P±) = 5(1 - x B )5 2 {p ± ). 



(101) 
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The fragmentation function 



D h/q (z, p hl _) = D (z, p hl _) 



1 



4z(2vr) 



dx + d x±e 



i(q x + + 



P h ± X J 



■) Tr <0\ip(x) a[a h ^(0)7~|0> 



=o, P h 



=o ' 
(102) 



is interpreted as the multiplicity of hadrons (of a certain type) found in the hadronization products (jet) 
of a quark with transverse momentum zero, hadrons with momentum fraction z = p^/k~ = p^/q^ 
and transverse momentum Ph±. Analogous to the case of the distribution function this can be seen by 
quantizing the quark fields on the lightcone, but now for x~ = 0. One then finds essentially a hadron 
number operator acting between two quark states. The sumrules discussed hereafter then confirm that 
D~ is a multiplicity distribution. If the original quark does have a transverse momentum, one should 
make the substitution p hl _ — ► p' hl _ = p hl _ — zk±, where p' h± is then the transverse momentum of 
the hadron with respect to the quark with momentum k±. Upon integration over p h± one finds 



D h/q{z) 



8tt 



dx e 



iq av 



Tr <0\ip(x)ala h ip(0)-f-\0> 



x-=x ± =o, p hA 



(103) 



This is the scaling fragmentation function. Both forms were first defined in |15] and also used in [16]. 
By integrating over z and using the techniques of appendix A we can derive the sum rule 



dz D h/q {z) = (n h ) q , (104) 
expressing the fact that D(z) is a multiplicity distribution. In [|l5| it was already proven that 



E 

hadrons 



dz zD(z) 



(105) 



which expresses momentum conservation. Finally we can also here consider the production of a quark 



from a quark. We substitute free quark operators for the aLa/j and find 

D q/q(z,Pq±) = 'H 1 ~ Z )^i.Pq±)- 



(106) 



The profile functions / and D 1 - do not have an interpretation as a distribution function, since 
they depend on interacting quark fields so no free field expansion can be applied. The anti-quark 
contributions can be derived completely analogously and lead to 



fq/H{x B ,p L ) 

D h /q(z, p h± ) 



1 



2(2vr) 3 
1 

4z(2vr) 



dx-d 2 x ± e^ q+x + P± X ^ Tr<P|V>(x)V>(0) 7 + |P> 
Ph±' x ± \ „—, , _ f 



dx d x±e 



i(q x + + - 



x+=0 

> <0|^(x) 7 -a[a h V(0)|0> 



(107) 



x~=0, p h 



(108) 



Inclusion of different flavours gives rise to a sum over contributing flavours and quark charge factors 



5.2 Semi-inclusive cross-sections 

Starting with the hadronic tensor derived in section 4, several cross sections for semi-inclusive lepton 
hadron scattering can be expressed in distribution and fragmentation functions. We will present a cross 
section for the process e + H —> e' + h + jet, i.e. not only a hadron is detected, but also the axis of the 
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jet it belongs to. Other cross sections will be determined from this one by integrating over measured 
variables. With the structure functions (|98|) and accounting for antiquarks and flavour we have 



^.(e+tf— e'+/i+jet) 

dx B dydzd 2 p hl _d 2 p ± 



8ira 2 ME 

Q 4 



E 



i=q,q 

+ 2(2-y)Vl-ycos 
1 



y + i-y 



x B fi/H(xB,P±)D h/i (z, p h ±-zp ± ) 



Ph±\ X B 

Q z 



fi/H(x B ,P±) 



\j; D h/i( z > Ph± ~ Z P±) ~ D h/i(z, P h± - zp ± )^ 



+ 2(2-y)y/l-ycos< 



\P±_\ XB 

Q z 



fi/ H {x B ,P±)D^ /i (z, p hl _ - zp ± ) 



[fi/H(x B ,P±) -x B f^ h (x B ,p ± )} D h/i (z, p h± -zp ± ) j. 



(109) 



Here <ph and 4>j are the azimuthal angles of hadron and jet with respect to the virtual photon direction. 
Note that p h± is defined with respect to the z-axis defined by the virtual photon, but the transverse 
momentum entering the fragmentation function p' h± = p h± — zp L is defined with respect to the jet 
axis. 

In our approach the transverse momentum of the jet immediately reflects the transverse momentum 
of the quarks within the target. To determine the jet cross section from this one has to first integrate 
over Ph±. A direct integration over p h± yields 



dx B dydzd 2 p L 



8na 2 ME 



E e f 



( 2 

y 



t=q,q 



+ 1-2/ 



2(2-y)^/l^coscf> j ^x 2 B f t j H (x B ,p 1 _)D h/i (z)\ (110) 



X B fi/H(x B ,P±)D h/i 



which gives for the jet cross section, 

da (e+H-,e'+h+3ct) ^C?ME 



dx B dyd 2 pj_ 



Q 4 



E< 

i=q,q 



( 1 



+ i-y 



X B fi/ H {x B ,PA_) 



2(2-y) v ^^cos^^x|/^(x B ,p ± )| (111) 



Note that we pick up a factor (nh) on both sides of the equation when integrating over z, which 
cancels. With this cross section we can calculate an expectation value for the azimuthal angle <^j e t, 



(COS </>jet) 



( 2\P±\ \ (2-yWT^y_ x B f^(x B ,p ± ) 
{ Q J 2(1 -y) + y 2 { f(x B ,p ± ) 



(112) 



where the structure factor given by the ratio of Xgf 1 - and / comes as an extension of the result 
previously derived by Cahn in [ 20 1 for the case of free quarks. Integrating over p ± in eq. ( |110|) we find 
the one-particle semi-inclusive cross section. The z dependence is simply given by a factor D(z), and 
the cosine integrates to zero due to the azimuthal symmetry of f ± (x B ,p ± ), 



da (e+H^e'+h) 8lTa 2 ME 



dx B dydz 



Q> 



E e i x Bfi/H(x B )D h/i (z). 



(113) 



i=q,q 
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If we compare this with eq. ([37]) we see that we have recovered the factorization and scaling of the 
semi-inclusive structure function 



H 2 {x B ,Q 2 ,z) = H 2 (x B ,z) = ^e 2 q x B f q / H (x B )D h/q (z) + f 7j/H {x B )D h/7j (z) , (114) 



and the semi inclusive Callan-Gross relation 

2x B H 1 (x B ,z) = H 2 (x B , z) 



(115) 



We recover the Callan-Gross relation because the corrections coming from the twist 3 diagrams are 
only entering in H3 and H4. Finally we will consider the cross section ( f35|) from which we started for 
I + H—*V + h + X. To this end we have to integrate over the jet transverse momentum. For the 
functions Hi occurring in eq. (35) we find 



Hi{x B ,z,Q 2 , p 2 h± ) 
H 2 (x B ,z,Q 2 , p 2 h± ) 

n 3 (x B ,z,Q 2 , p 2 h± ) 



2 E e * J d 2 P±fi/H(x B ,p ± )D h/i (z, p h± -zp ± ), 

i=q,q 

B H e< i J d2 P±fi/H(x B ,P±)D h/i (z, P h± -zp ± ), 



l=q,q 

2x B 

Z 



i=q,q 



E / d 2 P± 



P± ■ Ph± 

C i/h\ 



fi/h(x B ,p ± )D i/h (z, p h± -zp ± ) 



+ z [fi/h(x B ,P±) - x B f^ h (x B ,P±)) D i/h (z, p h± -zp ± ) 
+ fi/H(x B ,P±) ^D^z, p hl _ - zp L ) - D h/i (z, p hl _ - zp ± )J I 



(116) 



To find the last expression we used 



/ d 2 p 1 _p ti L f{x B ,p ± )D(z, p h ±-zp ± ) = f d 2 p ± (p ± - p h± )f(x B ,p ± )D(z, p h± -zp ± ). (117) 

We conclude that these structure functions do not factorize without extra assumptions on their trans- 
verse momentum dependence. However we do read off an extension of the Callan-Gross relation, 



H 2 (x B ,z, p 2 h± ) =2x B Hi(x B ,z, pl ± ). 



(118) 



6 Conclusions 

In this paper we have analyzed the structure functions in semi-inclusive deep-inelastic lepton-hadron 
scattering. We have considered all four structure functions, two of which (Wlt and Wtt) can only be 
measured by considering explicitly the momentum component of the produced hadron perpendicular 
to the virtual photon momentum, or by considering explicitly the perpendicular momentum of the 
produced jet. In all cases we have only considered the contribution of one (forward or current) jet 
being produced and particles therein. In the analysis we have included the twist two and twist three 
contributions. 

The twist two pieces in the structure functions can be expressed in terms of the well-known quark 
distribution and fragmentation functions which can be considered as specific projections of quark-quark 
correlation functions (profile functions). The twist three pieces, which constitute the main contributions 
in Wlt and Wtt involve two new profile functions that appear as projections of quark-quark correlation 
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functions. Up to order 1/Q, however, one must also include quark-quark-gluon correlation functions in 
order to obtain a gauge invariant result. By virtue of the QCD equations of motion the contributing 
pieces do not introduce new profile functions. The inclusion of the 1/Q contributions does not affect 
the evolution of the twist two profile functions (the quark distribution and fragmentation functions). 
We have not considered the evolution of the new profile functions, nor the 0{a s /Q) contributions in 
the process. An extension of the analysis to 0{Q~ 2 ) along the same lines as discussed here seems to 
be very hard to us. In that situation not only several new profile functions will appear but one has 
also lost the separation between kinematic quantities appearing in the distribution region from those 
appearing in the fragmentation region. 

Results in the on-shell parton model such as those previously derived by Cahn [^0| follow immediately 
as a special case by defining the profile functions with respect to a free quark target, i.e. the amplitudes 
in ( p8[) are then A\ = At = 4vr 4 5 4 (P - p) and B = 8ir 4 m q 5 4 (P - p). In this case the azimuthal 
asymmetry (costp) is of purely kinematical origin. 

The analysis in terms of quark correlation functions is of interest as it provides a convenient method 
to relate cross sections of hard processes directly to matrix elements of (nonlocal) quark and gluon 
operators. The latter could in principle be calculated when we knew how to solve QCD, or they can 
be calculated in a model. Although some data for semi-inclusive muon scattering exist [|2l|] we will 
not discuss them here. Neither have we presented any model as this would require extra assumptions. 
In the presented analysis a minimal amount of assumptions was made, namely those needed for the 
factorization of the hard scattering part. We also want to point out that one can do exactly the same 
analysis for Drell-Yan processes and e + e~-annihilation, thereby testing the validity of this approach and 
extending the experimental possibilities of measuring the profile functions. 

We would like to thank J.W. Bos, M. Dekker, A. W. Schreiber and R.D. Tangerman for stimulating 
discussions. This work is supported by the foundation for Fundamental Research of Matter (FOM) and 
the National Organization for Scientific Research (NWO). 



A Intermediate states and completeness 

The cross sections we are considering, describe the detection of one hadron in correlation with a 
scattered electron. In an experiment often many hadrons of the same type are produced in one event. 
This means that we have to pay special attention not only when we are integrating over the hadronic 
variables when we want to find an inclusive result, but also when we isolate an hadronic operator. In 
this appendix we discuss how to treat these situations, based on an example situation with only one 
type of scalar particles. Extensions to flavour and spin are obvious. An n-particle state is denoted 
as \pi, . . . ,p n > and we abbreviate the phase space integrations as dp = ^ 2 n) A 2E • ~'~' ie orc ' er °f tne 



particles appearing in the states has no meaning. Therefore the identity in this space looks like [g2| 



E 2 "' ( 119 ) 



n=0 

where 



X n = — ( dpi---dp n a\pi)---a^(p n )\0><0\a(p 1 )---a(p n ) (120) 

Th. J 



and 



|0><0|. (121) 



We first consider the situation with a multi-particle state where we integrate over all particles except 
one, the situation as we encounter it in our expressions. We denote with Px this complete system 
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minus one particle, 

d*P* 
(2ir) 3 2E x 1 

+ 2~! / dpidp2 \ pi ' P2 > Ph><PuP2 ' Ph \ + ■■■ 
= a} h la h = a[a h (122) 



/d ± f 
(27r) 3 2£ x \ Px ' Ph><Px ' Ph \ = \Ph><Ph\+J dp 1 \p 1 ,p h ><p 1 ,p h \ 



thus we can replace the sum minus one particle by the hadronic number operator. Then we want to 
integrate over p^ to reconstruct an inclusive formula. 



J (2Tr) 3 2E h (2tt) 3 2E x h>< = J dph \ ph><Ph \ + J dp hdpi\Pi,Ph><Pi,Ph\ 

^ J dp h dp 1 dp 2 \pi,P2,Ph><Pi,P2,Ph\ + • • • 



+ 2! 



(123) 



This is the identity operator, but with every term separately multiplied with the number of particles. 
The summation over all states then yields just the average number of particles h produced in the 
processes considered, 

/(fetfe'^^^?* 1 "' (124) 

Evaluated between particle states and reexpressed as cross sections this relation reads 

dn '' dE '- iahd wE C = Mn " E < :>) i£k- (125) 

B Structure of the correlation functions 

In this appendix we study the projections of the correlation functions defined as 

f A (P;x B , P± ) = -^Jdx'd 2 x ± e^ +x - + P^ x ^<P\lp(x)T A m\P> 



x+=0 



4zD B ( Ph ;z,k ± ) = J dx + d 2 x ± e^~ x+ - k ^ x ^ Tr <Q\^{x)a\a h ^{Q) r B |0> 



(126) 
(127) 

x~=0 



(with q + = —x B P + and q = p,/z). Using parity and time reversal invariance as well as hermiticity 
one proofs that the most general structure of the matrix elements is 

<P\^(x)^(0)\P> = J e { P- x [Ai(P,p)f + A 2 (P,p)P + B(P,p)]ji , (128) 

<0|^(x)4a ft ^(0)|0> = J-^e- ikx [C 1 (p h ,k)fi + C 2 (p h ,k) 1 f> h +V(p h ,k)] ij . (129) 

At this point one immediately sees that for unpolarized distributions f a and Dp are zero except for 
T A ,B = 7 M or T AyB = 1. 

Consider as the first case T A = 7^ for the /-profile functions. One obtains (choosing a frame where 

P±=0) 



f^(P;x B ,p ± ) = J dx~ d 2 x ± e i{q+x - +P ^ <P\^(x) 1 + m\P> 



x+=0 
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2(2vr) 3 
1 

W) 



dp 



dp 



2tt 

2 



Tr [ 1 + {Aij> + A 2 P + B)} 
A 1 (P,p) + —A 2 (P, P )} 

x B ) 



= r(x B ,P±), (130) 



fh-\P;x B ,p ± ) = 



J dx~ d 2 x ± e^ q+x +PrXi.) < p\ip( yX )~ / -ip(0)\p > 

{{p 2 + vl) 



x+=0 



2(2vr) 3 J 2vr 

M 2 
2(2vr)4(g+)2 



dp 1 



M 2 



-Ai(P,p) + x b A 2 {P,p) 



M 2 x\ 

2^)2 / &b,P±), 



(131) 



/^l(P;x s ,p ± ) 



2(2vr) 
1 



^ [ dx~ d 2 x ± e i{q+x +P± X ^ <P\4>(x)j ± 4>(0)\P> 

inr J 

dp 



x+=0 



2(2vr) 3 
P± 



2vr 



Tr [tj-^ + ^ + B)] 



(2tt)V 



fdp 2 A 1 (P, P )^-^f ± ^ B ,P±)- 



(132) 



Here / + ,/~ and are the profile functions. We assume that the functions A\ and B have no 
singularities, so that the profile functions are all of order 1. In the frame A = Q we also see at which 
twist these functions enter the calculations. When / + , i.e. the normal parton distribution function, 
enters at 0(1), then f 1 - enters at 0(Q~ r ) and / _ at 0(Q~ 2 ), as we expected from the discussion in 
section 3.1. It is straightforward to express as 



fW(P;x B ,p ± ) 



Q 



1 

Q 2 



(133) 



The function /W introduces another profile function e(x B ,p_ L ), 

fW(P;x B ,p ± ) = ^J— j dx-d 2 x ± e^ +x - + P^ <P\^(x)^(0)\P> 



2(2tt) 2 
1 

2(2^)3 
M 



x+=0 



dp 
~27 



Tr [Aij> + A 2 P + B] 



(2vr)V 



dp 



2 B(P,p) 
M 



Mx, 



e(x B ,p ± ), 



thus leading to the result 

f [1] (P;x B , P± ) = q ^ 2 -^e(x B , P± ). 



(134) 



(135) 



Although the profile function e appears to enter at twist 3, it will only give contributions at 0(Q~ 2 ) 
to the hadronic tensor. This is because e is always connected with the analogous Z)M function, which 
is also twist 3. 

In the case of a quark target we can explicitly calculate Ai and B and find 



f + (x BlPl _) = f (x B , Pl _) = f ± (x B ,p ± ) = e(x,p ± ) = 5(1 - x B ) 5 2 (p ± ). 



(136) 
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The analysis of the D-profile functions can be performed along similar lines. First again consider 
Ta = 7 M - In this case we want to consider the profile functions in the rest system of the produced 
hadron h. For that purpose we perform the following Lorentz transformation, 



z Ph 
[k~,k + ,k ± ] 

[D-,D + ,D ± ] 



[D-',D + ',D' ± ] 



2(Pfc) 



Ph 



ID', 



Ph 



2(Ph' 2 



Ph 



(137) 
(138) 

(139) 



where m 



hi 



m; 



+ Ph±- The profile functions in the primed system ( p' h± = 0±) are then given by 



(note that the hadron operators now have as argument p' h ) 

4zD^'(p' h ;z,k' ± ) = — 1— [ dx + d 2 x ± e^~ x+ - k '^ x ^ Tr <0|^(x)4o ft ^(0)7 _ |0 

(2ir) J 



> 



x~=0 



4zD^'(p' h ;z,k' ± ) 



(2vr) 
2 

1 



l — A Jdk + Tr [ 7 -( Cl f + C 2 A + n 



dk 2 



dx + d' 



k~=q 



k =q 



= 4zD-(z,-zk' ± )(U0) 



;(q-x+-k' ± -x±) 



Tr <0|V(x)4a fe V(0)7 + |0> 



(2vr) 3 

l — A jdk + Tr [ 7 +( Cl f + C 2 A + n 



x~=0 



(27T) 



fc =(jr 



(2vr) 4 (g-) 2 7 



1/;J ^)c 1 (^,fc') + -C 2 (^,A;') 



k =q 



4zD^'(p' h ;z,k' ± ) 



^3 J dx+d 2 x ± e i{q ~ x+ - k '^ Tr <0|^(x) 4o ft ^(0) 7 ±|0 

7^)4 / Tr [7x(Cir + c 2 A + P)]| fe _ =? - 

2k', 



> 



r k 

J dk 2 Ci{p' h ,k') =4^D ± (z,-zk' ± ) 



(2tt)V 

From this we can immediately find the values in the original frame, 

] (p h ;z,k ± ) = D~(z,p h± -zk±), 



D^ + \p h ;z,k ± ) = 



D^(p h ;z,k ± ) = 



2^f D+{Z ^ ~ ^ + 2z^FY 2 D ' {Z ' Ph± ~ Zk±) 

+ Ph± -%: )2 Ph±/z) vH Z ,Ph,-^) 

k± - p h ±/z 



zq~ 



D ± (z,p h± -zk ± ) + ^D (z, Ph± -zk_ 
zq 



In terms of T, q and the transverse momenta one has 

2{kl-plJz) 



D^(p h ;z,k ± ) = 



Q 



(f" + q^ D~ + 



zQ 



zQ 



1 

Q 2 



(141) 

x"=0 

(142) 
(143) 

(144) 
(145) 

(146) 
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The scalar function DW introduces the additional profile function E(z,p h± — zk±), 



4zD^(p h ;z,k ± ) = —^-3 f dx + d 2 x ± e i ( q * + +k ± -x ± ) Jr <0 |^( x ) a t aft ^(o)|0> 




dk+ Tr [Ctf + Ctfh + V] 



k =q 



2l7lh f ji2 V VPh,k) 4m h 
(2ir 4 <? J m h q 



(147) 



leading to the result 



D^( Ph ;z,k ± ) 



<£_ 2m h 
Q Q 



E (z,Ph± - zk±). 



(148) 



We find for the quark^quark profile functions 



D~(z, k ± ) = D + (z, k ± ) = D ± (z, k ± ) = E(z, k ± ) = 5(1 - z) 5 2 (k ± ). 



(149) 
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Figure 1: The general hadronic tensor (a) and the hadronic tensor with quark currents (b). 



Figure 2: The Born term in the expansion of the hadronic tensor. 



Figure 3: One gluon contributions in the expansion of the hadronic tensor. Only quark contributions 
are shown. 



Figure 4: Definition of scattering plane and momenta in IH — > I'hX. 



Figure 5: Gluon contributions giving rise to logarithms. 



Figure 6: Rules for the calculation of the diagrams. 



Figure 7: A general multi-parton contribution. 
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